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It is demonstrated, that chirality violating condensates in massless QCD arise from zero mode 
solutions of Dirac equations in arbitrary gluon fields. Basing of this idea, the model is suggested, 
which allows one to calculate quark condensate magnetic susceptibilities in the external constant 
electromagnetic field. 
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Consider QCD action in Euclidean space-time 



ipf{ij^'Vi_, + imf)ipf 



f '- 



where G^^^ is gluon field tensor, the sum is over quark flavours. 
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and is the gluon field. Pay attention, that in Euclidean formulation of QCD t/j is replaced by (The review 



of Euclidean formulation of QCD and instantons is given in see especially ■) The Dirac equation for massless 
quark in Euclidean space time has the form: 



where tpnix) and A„ are the eigenfunctions and eigenvalues of the Dirac operator —V 
fields operators into the left and right ones 

^+ = ^i^(l+75)+^'A^(l-75), 

where 

Then for nonzero A„ the Lagrangian and the action reduces to the sum of two terms 



(3) 



-i7^V^. Expand the quark 



L = - 



d^x 
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completely symmetric under interchange L < — > R. Therefore the solutions of the equations for left and right quark 
fields are also the same - the states, constracted from left and right quarks are completely symmetrical. This conclusion 
was obtained for fixed gluon field. It is evident, that the averaging over the gluon fields does not change it. Quite 
different situation arises in case Ao = 0. The contribution of this term to the Lagrangian: 

AL^y d^x[^j+ + 4'^]Vt^o (7) 

is equal to zero and no conclusion can be done about the symmetry of states build from left and right quark fields. 
One of the consequences from the said above is that all chirality violating vacuum condensates in QCD arise from 
zero mode solutions of Dirac equations 

These general arguments are supported by the well known facts: 

1. The general representation of the trace of quark propagator S{x) is expressed through the spectral function p{X) 
as a function of eigenvalues A (Kallen-Lehmann representation): 

TrS{x^) J dXpi\)A{x^,X) (8) 

At = A(x^, A) reduces to S{X) and we have (in Minkowski space-time): 

p(0) = -7r(0|?^(0)^(0) |0). (9) 

(The Banks-Casher relation ^). 

2. The zero- mode solution of ^ for massless quark in the instanton field is the right wave function - iPr{x) = 
(1 — l5)''P{x) and in the field of anti-instanton is the left one - iPl{x) = (1 -f 75)V'(a;) 4, ^. 

Basing on the statements, presented above, let us formulate the model for calculation of chirality violating vacuum 
condensates in QCD. Suppose, that vacuum expectation value (v.e.v.) of the chirality violating operator Oc.v. is 
proportional to matrix element iPqOc.v.iPo, where ipo is the zero-mode solution of Eq.Q in Euclidean space-time: 

(0 1 i^Oc.v.iJ 1 0) - Vf[a...V'o. (10) 



Ipo depends on x, on the position of the center of the solution Xc, as well as on its size p: tpo — ipoix — Xc, p)- Eq. fTT 
must be integrated over Xc, what is equivalent to integration over x ~ x^. (In what follows the notation x will be used 
for X — X(..) We assume, that p =Const and find its value from comparison with the known v.e.v. 's. Finally, introduce 
in pop the coefficient of proportionality n. So, our assumption has the form: 



(0 I ^(0)Oe...V'(0) \0)=-nJ d^x^+{x,p)Oc.v.Mx,p) (11) 

Our model is similar to delute instanton gas model jo], where Xc is the position of instanton center. Unlike the latter, 
where the instanton density has dimension 4, n has dimension 3 and may be interpreted as the density of zero- modes 



centers in 3-dimcnsion space. Note, that the left-hand side of pT|) is written in the Minkowski space-time, while the 
right-hand side in Euclidean ones. (The sign minus is put in order to have n positive.) For x and p-dependens of 
ipo{x,p) we take the form of the zero-mode solution in the field of instanton in SU{2) colour group: 

M^:P) = 1(^-15)1 j^T^^^^Xo, (12) 

where xo is the spin-colour isospin (| T |= 1/2) wave function, corresponding to the total spin I + T = J equal to 
zero, J = 0. tpoix,p) is normalized to 1: 

J d'^xi;+{x,p)i>{x,p) = 1 (13) 

Consider first the quark condensate (0 | | 0), the most important chirality violating v.e.v.. determining the 
values of baryon masses . (Here q = u,d are the fields of u, d-quarks). In this case Oct,. = 1 and in accord with 
([n]),(I13l) we have 

n = -(0 I I 0) = (1.65 ± 0.15) x 10"^ Q^yS ^^^^ ^ Q^y^ j^qj f^^^-^ 

(The integration over SU{2) subgroup position in SU{3) colour group as well as anti-instanton contribution are 
included in the definition of n.) The anomous dimension of quark condensate is equal to 4/9.) According to (|14|) n 
has the same anomalous dimension. The size p of the zero-mode wave function can be found by calculation in the 
framework of our model of the v.e.v. 

- 3(0 I ^T^.^G^i, I 0) EE mliO \qq\0). (15) 



The parameter ?71q is equal to 



11|: ?71q = 0.8 GeV^ at 1 GeV. The tuq anomalous dimension is equal to -14/27. 



Working in the SU{2) colour group, substitute A" by T°(a = 1, 2, 3) and take for C^^jj the instanton field 

4 p2 

G%{x,p) = -??ap^ , „ (16) 

n n 

where the parameter rja^iy were defined by 't Hooft [IJI (see also [0|). The substitution of and into pT|) gives 
after simple algebra 

= mln. (17) 

2 

Therefore, 

p = = 0.79 GeV"^ = 0.156fm(at 1 GeV). (18) 

V 2too 

We are now in a position to calculate less well known quantities - the magnetic susceptibilities of quark condensate, 
induced by external constant electromagnetic field. 

The dimension 3 quark condensate magnetic susceptibility is defined by 



13| 



{0\qa^,q\0)F ^egx{0\qq\0)F^,, q = u,d, (19) 



where quarks arc considered as moving in external constant weak electromagnetic field i^^^ and eg is the charge of 
quark q in imits of proton charge (the proton charge e is included in the definition of -F^^). The left-hand side of 
(frojl violates chirality, so it is convenient to separate explicitly the factor (0 | | 0) in the right-hand side. It was 
demonstrated in 1^ that (0 | q<Jnu<l | 0)_f is proportional to the charge eg of the quark q. A universal constant x is 
called the quark condensate magnetic susceptibility. 

Let us determine the value of x in our approach. For this goal it is necessary to consider Eq[3] in the presence of 
external constant electromagnetic field F^j^i, and to find the first order in F^^ correction to zero mode solution (jl2p . 
This can be easily done by representing t/i as 

^{x,p)=ipo{x,p)+i)i{x,p), (20) 

where -0o is given by ([T^ and ■01 represents the proportional to F^^, correction. Substitute (PP]) in EqO added by the 
term of interaction with electromagnetic field, neglect ij^i in this term and solve the remaining equation for ipi(x, p)). 
The result is: 



1 f lx'^\ 

'^i{x,p) = jQegrjai^i.craFfj.^x'^ \1 + -— \ipo{x, p), 



(21) 



where CTq are Pauli matrices. The matrix element '4'~^(T^i,ip appears to be equal 



1 / Ix^' 

ib+a^,ij = ^-egF^,i4x^l + -^ ]^o- (22) 



(The properties of iiat^iy symbols [1^,01 were exploited.) The v.e.v. in the Minkowski space-time is given by: 

(0 I ^a,,^ I 0)^ = e,F,^n^ J d^xx' (l + J^^^^ (23) 

(The normalization condition (|13|1 for ^o{x, p) was used.) It is convenient to express n through quark condensate by 
([H]) . use the notation = r'^, where r is the radius- vector in 4-dimensional space. Then according to we have: 

X^V/*V(1 + 1I;)^^ (24) 



The integral is quadratically divergent at large r. So, the cut-ofF R is introduced. Its value can be estimated 
in following way. The volume occupied by one zero- mode in 3-dimensional space is approximately equal to 1/n (the 
volume of the Wigner-Seitz cell). So, for cut-off radius square in four-dimensions we put 

where the factor 4/3 corresponds to transition from 3 to 4 dimensions. The calculation of the integral (pi)) at the 
values of parameters p (fT8|) and R^ ([25|) gives 

T~2 



X = -3.52GeV~^ (26) 



The quark condensate magnetic susceptibility was previously calculated by QCD sum rule method [ij] 
expressed through the masses and coupling constants of mesonic resonances. The recent results are: 



and 



x(lGeV) = -3.15 ± O.SGeV"^ [16]; x(l GeV) = -2.85 ± 0.5GeV"^ [17] 

n 

(The earlier results, obtained by the same method, were: x(0-5 GeV) =—5.7 GeV ^ jl4J] and x(l 



(27) 



-4.4±0.4 GeV-2 



n 



GeV) 



171 ]. but not in 



15|. The anomalous dimension of x is equal to -16/27. It was accounted in [l4] 
the presented above calculation. (In some of these papers, the as-corrections and continuum contribution, were also 
accounted.) One can believe, that the value (QH) refer to 1 GeV, because the value of quark condensate refer 



to this scale and also because the scale 1 GeV is a typical scale, where, on the one hand, the zero-modes and quark 
condensates are quite important (see, e.g. 10]) and, on the other, the instanton gas model is valid jfi]. Since the 
integral is quadratically divergent it is hard to estimate the accuracy of ([26|) . I guess, that it is not worse, than 
30-50%. In the limit of this error the result p6p is in an agreement with those found in phenomenological approaches. 



Turn now to quark condensate magnetic susceptibilities of dimension 5, n and ^ defined in Ref. 

g{Q\q\rGl,q\Q)F^e,KF^,{0\qq\Q), 



(28) 



*5e^,p,(0 ] q^z-V'G^p^q \ {))f = e,^F^,{0 \ qq \ 0) 



(29) 



Perform first the calculation of k. In this case the expression of ipi{x,p) ()2ip must be multiplyed by the additional 
factor: \t^G^ where G^^^ is given by and the indices n, v in (PT|) are changed to A, <j. In the further calculation 



it will be taken into account, that xo hi (|12[) corresponds to total spin-colour isospin J = and consequently 

TTXo = -3X0, ctV^xo = -'^'^''Xo. (30) 



In the relation 



Vbfj.uVbXcr = S^xSi^rr — Sf^^SuX + £^iv\o 



(31) 



the last term drops out after summation of zero-modes from instanton and anti-instanton configuration. The final 
result for k is: 



1 



■"^du- — \ 1 -I- -M 



{u + If 



w .^ 13 1 11 

ln(z + l)- + + 



1 1 



6 z + l 2 (z 4-1)2 3(^ + 1)3 



where z = R / p = 12.7. Numerically, we have: 



K = -0.26 



(32) 



(33) 



The calculation of ^ is very similar to those of k, and the result is 



^ = 2k = -0.52 



(34) 



The values of k and ^ only logarithmically depend on the cut-off. But unfortunately the logarithm in ([32)) is not 
very large and its main part is compensated by the term —13/6, appearing in p2p . So, the accuracy of (|33p . (|34p 
can be estimated as about 30%. The phenomenological determination of 5-dimensional quark condensate magnetic 

:in 



susceptibilities was performed by Kogan and Wyler 



dimensions were accounted. The results of 



18l | along the same lines, as it was done in [1J1,[15(. No anomalous 



181 are: 



/t== -0.34 ±0.1, ^ = -0.74 ±0.2 (35) 

As can be seen, they are in a good agreement with ([33)) . ((34|) . The 5-dimensional quark condensate magnetic suscep- 
tibilities play a remarkable role in determination of A-hyperon magnetic moment 19[ . 

I conclude. It was argued, that chiral symmetry violation in QCD arises due to zero-mode solution of Dirac equation 
for massless quark in arbitrary gluon field. The model is proposed similar to dclutc instanton gas model, in which the 
zero-mode solution is the same as in the field of instanton. The parameters of the model: the density of zero-modes 
and their size are determined from the values of quark and quark-gluon condensates. In the framework of this model 
the values of quark condensates magnetic susceptibilities of dimensions 3 and 5 were calculated in agreement with 
ones found by phenomenological methods. The success of the model supports the basic idea of our approach and 
shows that it can be used in determination of other chirality violating condensates. 

This work is supported by RFBR grant 09-02-00732. I acknowledge the support of the European Community- 
Research Infrastructure Integrating Activity "Study of Strongly Interacting Matter" under the Seventh Framework 
Program of EU. 



[1] Instantons in Gauge Theories, ed. by M.Shifman, World Sci., 1994. 

[2] A.I.Vainshtein, V.I.Zakharov, V.A.Novikov and M.A.Shifman, Usp. Fiz. Nauk 136, 553 (1982) 
(Sov. Phys. Uspekhi 25, 195 (1982)), Reproduced in [1], p.468. 

[3] T.Banks and A.Casher, Nucl. Phys. B169, 103 (1980). 

[4] L.S.Brown, R.Carlitz and C.Lee, Phys. Rev. D16, 417 (1977). 

[5] R.Jackiw and C.Rebbi, Phys. Rev. D16, 1052 (1977). 

[6] T.Schafer and E.V.Shuryak, Rev. Mod. Phys. 70, 323 (1998). 

[7] B.L.Ioffe, Nucl. Phys. B188, 317 (1981), Errata 192, 591 (1982). 

[8] B.L.Ioffe, Zs. Phys. CIS, 67 (1983). 

[9] B.L.Ioffe, ArXiv: 0810.4234, to be published in Phys. At. Nucl. 72, No. 7 (2009). 

[10] B.L.Ioffe, Prog. Part. Nucl. Phys. 56, 232 (2006). 

[11] V.M.Belyaev and B.L.Ioffe, Sov. Phys. JETP 56, 493 (1982). 



7 

[12] G. 't Hooft, Phys. Rev. D16, 3432 (1976). 

[13] B.L.Ioffe and A.V.Smilga, Nucl. Phys. B232, 109 (1984). 

[14] V.M.Belyaev and I.I.Kogan, Yad. Fiz. 40, 1035 (1984). 

[15] I.I.Balitsky, A.V.Kolesnichenko and A.V.Yung, Sov. J. Nucl. Phys. 41, 138 (1985). 

[16] P.BaU, V.M.Braun and N.Kivel, Nucl. Phys. B649, 263 (2003). 

[17] J.Rohrwild, JHEP 0709:073 (2007). 

[18] I.I.Kogan and D.Wyler, Phys. Lett. B274, 100 (1992). 

[19] J.Pasupathy, J.P.Singh, S.L.Wilson and C.B.Chiu, Phys. Rev. D36, 1442 (1987). 
The end of submission. 



